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1 Introduction 



In HJI we have proposed a quark model of current matrix-elements in the heavy quark limit 
which has been shown to present two important features : it is covariant and it presents the 
full scaling properties of heavy quark symmetry (HQS), demonstrated by Isgur and Wise 
for QCD 0. The model is based on an old formulation of relativistic multiparticle states, 
i.e. a representation of the full Poincare group, proposed by Bakamjian and Thomas 
and reexpressed by Osborn Q . We have given a simpler and concise description of Poincare 
generators and wave functions in momentum space Q. Since some time this formulation, 
or variants of it, have been used by a large variety of people to formulate quark models of 
form factors, especially on the null-plane ||. One must stress that by boosting the states to 
arbitrarily large momenta, one recovers the models directly formulated on the null-plane. 
In particular, one gets null-plane wave functions with the correct kinematical limits, unlike 
other approaches. 

It must be emphasized however that in this formulation, for finite quark masses, 
the currents based on the free quark current operators are not covariant. What we have 
found is that their heavy quark mass limit are covariant. One important consequence of 
this covariance property of the limit is that all the models we have referred to will have 
essentially the same limit, whichever the frame in which they have been formulated when 
rriQ < oo. Then, the only remaining source of variation consists in the choice of the 
mass operator or of the set of wave functions at rest, and everything we are to say, being 
completely independent of such choices, should hold for the rrtq — > oo limit of such models. 

The implementation of covariance and heavy quark symmetry obtained in this ap- 
proach must be fully appreciated. Indeed they are not built in, neither enforced by hand. 
It represents an important progress with respect to older popular models of heavy-to-heavy 
semi-leptonic form factors which, on the one hand, make implicit or explicit reference to 
particular frames to calculate form factors, and, which, on the other hand, finally renounce 
to predict form factors except for a privileged q 2 value, and inspire for instance from VMD 
to extrapolate form factors. Then either heavy quark symmetry is not satisfied or it is 
enforced by hand. In the present approach, on the contrary, everything results from a 
systematic calculation from the eigenstates of a mass operator, and a current operator. It 
must be also underlined that the approach, while leading to such nontrivial properties, does 
not spoil the simplicity which makes quark models so attractive. It preserves principles of 
quark models like a fixed number of constituents, a three-dimensional description through 
ordinary wave-functions, a free quark current operator . On the other hand, it is also 
a progress with respect to our own older models which, while presenting heavy quark 
symmetry, were only approximately covariant (i.e. the Isgur- Wise scaling function £ was 
depending on the chosen frame). 

In the present paper, we pursue the investigation of the model presented in QTJ| 
as regards general properties, not dependent on particular choices of interaction or wave 
functions at rest : the mass operator will remain totally arbitrary, except of course for 
rotation and parity invariance. As another important advantage of the new approach, 
we want now to show that it may help a lot in understanding more physically, in the 
context of bound state physics, sum rules already formulated in a field-theoretical context. 
Let us stress that exact saturation of sum rules is a rather specific feature of this type 
of relativistic models. Sum rules strongly rely on completeness relations. But what is 
precisely needed is completeness relations for wave functions of states in motion. Such 
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relations are automatically provided by the construction of a unitary representation of the 
Poincare group. On the other hand, many phenomenological models use wave functions in 
motion which are not shown to be orthonormal, e.g. to be eigenstates of an Hamiltonian 
or even, as recalled above, use form factors not calculated through wave functions except 
for a normalisation at some q 2 . Then, they have no reason to satisfy sum rules. 

In addition to demonstrating that Bjorken-Isgur-Wise (BIW) sum rule is exactly 
valid in our approach, we shall also decompose the contributions of the various states 
to the "physical side" of the sum rule into different parts with clear physical meaning, 
clarifiying in particular the nature of relativistic corrections and the origin of the various 
bounds on p 2 . 

2 Direct demonstration of duality and saturation of 
the Bjorken-Isgur-Wise sum rule 

The sum rule under consideration writes, for mesons, in the transparent formulation of 
Isgur and Wise [§] 

P 2 -E(hyi) 2 + 2T 3 fc /2 (l) 2 ) = i (i) 



where p 2 is the slope of the elastic ground state Isgur- Wise scaling function, while the Ti/ 2 , 
T3/2 are the scaling functions of the transition between the ground state and the P wave 
(L = 1) states with possible radial excitation, k — 1 for the lowest P state, k = 2, . for its 
radial excitations. In the following, the n's are always considered at the zero- recoil point 
w — 1. Therefore, to reduce notation, we denote the corresponding values by T1/2, t 3 / 2 . 

To demonstrate the sum rule, we shall not use the original Bjorken method start- 
ing from commutators of field theory M, since our model does not contain a priori such 
commutation relations. Rather, we shall calculate directly the sum of squares of transition 
matrix elements to any final state and find that it has a very simple value, independent 
of the dynamics. In the present approach, states have a fixed number of constituents N 
(no pair creation). Moreover, for a confining interaction, they will be stable bound states ; 
this implies that saturation of the sum rule holds within the resonance-dominance and 
narrow-resonance approximations. 

For the particular purpose of this article, it is found easier to treat on equal footing 
the various states with different spins and which may be mesons or baryons or states with 
any number of quarks ; we then abandon the manifestly covariant formalism with Dirac 
spinors, and return to the initial bidimensional-spinor formalism. In this formalism, current 
density matrix-elements write, without assuming m 1; — ► 00 : 

<p'\f\p>= Jud 3 Pt £ r s[Sl (p'-^m) 

i sisis[ 

< PAirwiS! > ^ slSl {p - zti, {pi})- (2) 

Index 1 denotes the active quark, while spectator quarks, from 2 to N, are denoted by the 
generic index i. We disregard color and flavor. The ip's are what we call the "relative" 
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wave functions. They are obtained by dropping 5(Hpi — P) in the total wave function with 
total momentum P (we normalize any state to S(P' — P) and 1 for any discrete index). 

Part of our demonstration will be worked out without passing to the heavy quark 
limit. Let us recall that at this stage the model is not covariant, but we do not need the 
covariance property. For the first and most general demonstration of the sum rule, we do 
not need either to make explicit the construction of the t/>'s from the internal wave functions 
at rest. This construction shall be explained in detail in section 3, where one really needs 
it. We need only to know that the ?/>'s may be chosen to form a complete orthonormal 
basis of the subspace with P fixed. Let us label it by a generic label n in addition to P, 
we mean that : 

E / T[ d3 Pi ^n's lSi (P ~ S Pi> {Pi})i>ns lSi {P - ^Pii {Pi}) = <W • (3) 

Then we have the closure relation : 

Etfc.^-EiS.GSHt^ • (4) 

n i i i 

It is of course convenient to choose the basis among the eigenstates of energy 
\J M% p + P 2 , and to choose for the label n the labelling of the corresponding internal wave 

functions at rest tp n (kx, {ki}), which will include, in addition to various internal spin and 
angular momentum labels, excitation numbers. 

Let us now consider the transition matrix-elements from one fixed state denoted as 
n = 0, P, to all possible states, i.e. with all possible n, with another common momentum 
P', induced by We choose for j^ 1 the elastic vector current, although we could choose 
any other current and even any Dirac matrix O, with mi 7^ m' 1 - indeed, it is known that 
whatever the choice, we would obtain the same Bjorken sum rule because of HQS. Our 
choice is made for transparency and definiteness only. 

With the closure relation eq.(^j) at hand, one can perform the following sum : 

h^(P, P')=E< n, PV|0, P > x (< n, P'\f\0, P >)* (5) 

n 

which we will call the "hadronic tensor" in analogy with inclusive leptoproduction or semi- 
leptonic decays. Let us stress that this sum should not be covariant in principle, even if 
the theory were covariant, because each intermediate state has a different energy at fixed 
P' . We find trivially from eqs. ([| and |]) : 

r N 

hr(p, P')= e j II d *p* ^ks+P' - E^» {pMos^p - te» 

(u si (p' - mru si (p - spo) x (u s[ (p' - mYu 7l (p - mr ■ (6) 

The important step already realized is to have reexpressed the hadronic tensor in terms of 
the initial state (n = 0) wave function only. In fact, the right-hand side of the equation is 
the average of a rather simple operator. Introducing the Fourier transform of the current 
density : 
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j(q)= d 3 xe-^ s j(x) 



(7) 



one can write : 



h^(P,P>) =< 0;Pj1/(0) 0;P> 



(8) 



with q = P' - P. 

Up to now all the written relations were valid in full generality. Let us now take 
advantage of the heavy quark mass limit m\ — > oo. The important point is that in this limit, 
the velocity of the active quark can be considered equal to the velocity of the corresponding 
hadron, v or v'. Then we can factor out of the integral the spinor factors, which now depend 
only on v and v' : 

h^{v^)= E [u^^UsM] [u s ^)Yu 7l {v) 



r N 

E / n d 3 Pi r 7 is s p - ^ te))^w(p - m, m) 



(9) 



A further simplification can be obtained if one averages as over the polarisation of the 
initial state ; let us call it A. Then it is easy to convince oneself that : 



N 1 
i=2 ZJ -T J- su x 



(10) 



First, this is easily seen to follow from rotation invariance if P = 0. Then, one can 
extend the result to arbitrary P using the following property of the model : in passing to 
wave functions in motion, the spectator spins s\ are simply rotated by Wigner rotations, 
which are cancelled by contraction, while for the active quark spin Si,Si, the Wigner 
rotation reduces to unity in the m\ — > oo limit. Then 



• (11) 



Whence our central relation : 



l. free quark 



(12) 



We have exact duality for the hadronic tensor, that is the general Bjorken sum rule 
in the form of Isgur and Wise. Let us recall that it holds with any Dirac space matrices 
instead of 7^, 7". 

Though it may appear trivial, it is worth stressing once more that such an exact 
duality can hold only because current matrix elements are actually calculated as matrix 
elements of a given operator between wave functions in motion that satisfy closure. If, on 
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the contrary, form factors are constructed according to phenomenological recipes, although 
one may have at start or t honor mal wave functions at rest, one cannot get duality for 
arbitrary q 2 or w = v-v'. Two properties of this heavy mass duality must also be underlined. 
First, the limit of h^ v is manifestly covariant, in contrast to h^ v at finite rriQ. Second, this 
limit is the sum of the heavy mass limit of each product of two matrix elements. This 
nonobvious exchange of limit and sum can be shown to hold within the model, using 
closure for the wave functions in the heavy mass limit. That entails that the sum of limits 
is convergent, which should not be necessarily true in general in field theory. 

From this duality relation, one deduces easily the sum rule (Q) by reexpressing the 
current matrix elements in the definition of h^ u {P,P') (eq.([5|)) in terms of the standard 
invariant scaling functions, . The deduction will be made by expanding /i o to second 
order around v, if = ; for simplicity we choose v and if collinear. Moreover, we take 
as initial state a 0~ meson. One finds, taking into account that our own normalisation of 
states is 5 3 (P — P'), different from the one of Isgur and Wise : 



At this point, we make use of the fact, demonstrated in the preceding letter, that the 
model indeed satisfies Isgur- Wise scaling, thanks in particular to the factorisation of the 
wave function at rest into spin and rotationally invariant space wave functions. This result 
extends to radial excitations, by dropping the 1 in eq. (|T3|), because of orthogonality of 
wave functions at rest. On the other hand, we have not yet presented the demonstration 
of covariance and scaling for the general case, which we postpone to another paper. We 
assume for the moment this general result. Manifestly covariant and scaling expressions 
will be presented in a forthcoming paper for the transitions to L — 1 (P wave states) 
Tofl . CT radial excitations yield 0({v' — if) 2 ) amplitudes, therefore a negligible 0({y' — v) A ) 



contribution to /loo- The other possible final states through second order in v, if in the sum 
rule are 1~ and the L = 1 + , 1 + , 2 + states, because j° is scalar under rotations. — > + 
vanishes identically for a vector current. — > 1~ and 2 + vanish since for j° they have the 
form v x if from covariance. Therefore, in addition to the elastic — > contribution, 
we have just the — > 1 + one, to express in terms of the scaling functions at zero recoil 
Tj. This is readily done : 

<l/2,e x *\f\0->*-T l/2 e x *-(if-v) (A = 0) (14) 

<3/2,e A *|j |0->^-v / 2r 3/2 £ A *-(t; , -t;) (A = 0) . (15) 

We have made an expansion only at the lowest, first order in v, if, since these expressions 
appear in squares (therefore we have made vq,v' q ~ 1). Then : 



£|<l/2,A|j°|0-> +^|<3/2,A|j |0-> = (|r 1/2 | 2 + 2|r 3/2 | 2 ) (v - v') 2 . (16) 

A A 

Finally, summing on the P-wave states, including excitations : 

h 00 (v,if) « 1 - p 2 (v- if) 2 + £ (\r k l/2 \ 2 + 2|r 3 fc /2 | 2 ) (v- iff . (17) 
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On the other hand, the right-hand side of eq. (|1^) is 



Ul eequark ?zl-^{v-v'Y (18) 

whence the desired sum rule eq. (|1|). 

For baryons, the difference will be that the scaling functions for the elastic transition 
are defined through the coefficient of 7^ instead of + v' 11 , whence : 

< 1/2VIV2+ >~ 1 - \p\{v-$? - \iy-$y (19) 

because : 

%7%~ ±-l(v-v>) 2 (20) 

o 

with our normalisations (u^ub = 1)- In this case, the sum rule does not contain the 1/4 
term, but this has no deep physical meaning : 

p 2 B = inelastic contributions . (21) 

One must note that one is dealing with the usual Bjorken sum rules, corresponding to 
the first term in the current commutators considered by Bjorken ("direct" contribution) . It 
has been shown by Bjorken that the direct and z-graph contributions satisfy two separate, 
independent sum rules J7|. In our model as it stands, the z-graph contribution is absent 
due to the absence of pair creation or annihilation, therefore the second sum rules are not 
present. But a slight modification of the model, including pair creation or annihilation by 
the current and extending correspondingly the space of states, would probably allow to get 
the commutation relations and to satisfy the second type of sum rules corresponding to 
the z-graph contribution. 



Analysis in terms of internal wave function matrix 
elements 



Now we pass to a more detailed calculation of the various contributions to the left-hand 
side of the sum rule, as expressed in terms of internal wave functions at rest (p, eigenstates 
of the mass operator. That is, we display the various effects due to the hadron center-of- 
mass motion as treated relativistically. Then, we get more physical insight. We use the 
expression (12) of the preceding paper 0, obtained by explicitation of the relative wave 
function in terms of yj's (this is before passing to the tjiq — > 00 limit) 



<P'\f\P>= / l[d 3 p 



N 



X 



i=2 



\ M' M 11 v 9^o 



x 



N 

EE^'({S» \D(R';\ llsl/ll u sllll Yus''MRi)s^ n^(< 7l ^k^({^}) • ( 22 ) 



i=2 
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The ki and k'^s are the internal momenta corresponding respectively to the initial and final 
hadrons. They are complicated functions of the spectator momenta pi, obtained through 



the equations (4) of the preceding paper |] . an d k® denote simply J mf j + p^ and 



mf + kf and M = y^kf, j runs from 1 to N. Rj are Wigner rotations defined in the 
same paper. Let us recall how these Wigner rotations are practically related to boosts. 
Let A be a Lorentz boost : 

\/p° A u s (p) = y^Ap)^ D ss , (R) u sl {Ap) (23) 
with the rotation R defined by : 

R = B$AB P (24) 

B a is the boost which applies (Vr, 0) on the four vector a M . The root factors in the 
formula are needed because of our normalization of states. From this formula, one can 
deduce easily D SS /(R) knowing A and p. 

Let us identify the origin of the various factors involved in eq. (|22|) , in addition 
to the wave functions <^'s and to the free quark current density. In front of ip', there 
is the product of square roots of the Jacobians corresponding to the change of variables 
pi) {Pi} ~^ {h}, P and p[,{pi} — > {ki},P' ; the presence of these factors ensures unitarity 
of the corresponding functional transformations. Between <p' and (p, the Wigner rotations 
of quark spins ensure the passage from the ordinary one-particle spins to the internal spins 
; note that the D matrices are unitary by themselves. 

In the limit uiq = wi\ or m' x — > oo, the expression simplifies considerably. For 
convenience, we choose v and v' along the same axis Oz and count them by the algebraic 
numbers v z and v' z . We maintain v ^ to check the requirements of covariance on matrix 
elements, as imposed by the expressions (|I3"D, (0), QI3D- First, the relation between the 
kiS and p^s becomes : 

fc? = v °p° i - v z p iz , k iz = v°p iz - v z p°, k iT = pi T , (25) 

where T denote the component perpendicular to Oz. One has an analogous relation for the 
k"s with v' instead of v. Then, the Jacobian factors take the very simple form : 



,-, y V°Pi - VzPizy/V^Pi ~ v' z Piz 
i Pi 

The Wigner rotations of the active quark tend to unity because its momentum becomes 
parallel to the direction of the boost, which becomes v or v'. Finally, the active quark 
current density tends to : 



u s [ 7° u Sl 



(27) 



and can be once more factored out of the integral. Whence the final expression in the limit 
m Q — > oo 



. N N 

< p'\f\p >= £ u s ^ Usi / n d 3 Pt n 



\J(Pi ■ v)(pi ■ v 1 ) 



i-2 i-2 Pi 
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N 



(2f 



i=2 



Let us now expand the various factors in this expression around v, v 1 — 0, at fixed pi, ki 
being a function of Pi and We note that at v — v 1 — 0, one gets simply the scalar product 
of the internal wave functions at rest (p, (p f . Therefore it is 1 for the elastic transition 0^0 
and for inelastic transitions. 

Then, to calculate the r» 's at zero recoil, we have once more to calculate inelastic 
amplitudes only through lowest, i.e. first order 0(v), or 0(if). Let us first calculate these 
inelastic contributions. We have the following lowest order expansion of the integral : 



<n|i°|0>= JU^PiH^iPi 




[a x p iT ) 



(29) 



Pi + rrii 

where scalar product on spin space is implied. The first term comes from the combination 
of the variation of the Jacobian factors, and the variation of the argument k of the wave 
function, the second one from Wigner rotations. Each has separately a factor v' z — v z , as 
required by covariance. Indeed, the covariance requires a factor (v' z — v z ) for the sum ; 
now, if we had dropped the spin and considered scalar quarks, the first term only would 
remain ; covariance then requires a factor v' z — v z on this factor separately. Note that in 
this approximation the active quark density u si r y°u si ~ 1 does not give a contribution. We 
can replace i d/dp iz by the more suggestive notation ; it is indeed the space coordinate 
operator. In addition, we henceforth denote matrix elements between ip n , ip n i as (n'||n). 
Then : 



< n\j°\0 >^ (v' z -v z 
If we particularize to mesons, i = 



E 



71 



p\izi + izip^ i (a; x p iT ) 



2 2 pj + m 

2 only ; we can drop the index i and write : 







(30) 



<r P \i°\ c t ~~>i 

^ 1 wave\J \ u wave 



t> - V, 



71 



p°iz + izp° i (<t x p T ). A 



2 p° + m 







(31) 



We now particularize to a initial state, and to n = 1 + . Then, one finds from the 
identification with eq. (O): 



1/2 



3/2 



E E 

1/2,3/2 k 



l + ,k 



p°iz + izp° i (a X pr) 



2 p° + m 



0" 



(32) 



It is clear that the operator on the right-hand side leads only to transitions to 1 + (it has 
L — 1). Therefore we can as well replace the summation on states by a summation on the 
whole Hilbert space, and using closure of the (p% we end with : 



E 

k 



r l/2 



+ 2 



T- 



3/2 



p z + zp 



+ 



x Pt) : 
p° + m 



(33) 
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We have taken into account that a average on S = states is 0, so that the two contribu- 
tions add in squares. Further simplification is obtained by using rotation invariance in the 
Wigner rotation contribution : 



E 

k 



1/2 



+ 2 



T 3/2 



0" 



' p°z + zp°" 



+ 



V 



6 (p° + m) 2 



0" 



(34) 



which is an average on space wave functions only. 

As to p 2 , we could simply refer to the preceding paper, eq. (29) of the published 
version]!]]. However, it is more instructive physically to recalculate p 2 along the same lines 
as we have just done for the r's, i.e. directly in the formalism with bidimensional spin and 
Wigner rotations. We have just to push the expansion in v and v' of the matrix element 
through second order. It happens that the result can be decomposed in a manner similar 
to the r's. One has three types of effects which make the matrix element depart from its 
zero recoil value, 1 : 

i) one from the Jacobians and the spatial wave function arguments variations 

ii) one from the Wigner rotations 

iii) one from the current the active quark, which was not present for P waves. 

noi Covariance requires the whole result to be of the form 1 — (p 2 /2)(v f — v) 2 . But 



in fact each effect gives a separate contribution, which is oc 



Indeed for iii) it 



is obvious, see eq. (j3~?|) below (it is the requirement of covariance for a one quark state). 
Then the effects i) and ii) must also combine to give a contribution of this form. Now, there 
are no crossed terms between them ; these should correspond to the product of factors of 
first order in the velocity : first order terms from the Wigner rotation effect ii) contain a 
spin operator a, which averages to on a 0~ state ; then, these first order terms cannot 
combine with one from i) (the latter effect does not generate spin operators), i) and ii) give 
therefore non-interfering additive contributions. Moreover each one must have the form 
{y 1 — v) 2 . For the first effect i), it is seen by dropping spin and considering scalar quarks : 
then i) would give the only contribution and therefore it must be oc (v' — v) 2 separately by 
covariance. It results that the same form must holds for ii). The contribution i) to p 2 is 
found to be, by a somewhat lengthy calculation : 



P. 



space 



p z + zp 



(35) 



The contribution ii) is 



P wigner 

Finally, from the active quark current 



1 p 



6 (p° + my 



(36) 



<pWi\f\PiSi >= 1 - i-itf - v" 2 



one has : 



(37) 



(38) 



whence : 
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p 



2 

P space 



+ P 



wigner 



+ 



p z + zp 



+ 77 



P 



6 (po 



m 




(39) 



BIW sum rule is obviously satisfied, considering eq. 

It is seen that "spatial" contributions to p 2 and r 2 's, i.e. the one obtained by 
dropping quark spin, cancel each other as it should since, for scalar quarks, the sum rule 
should write : 



P 



E 



T 







(40) 



since in this fictitious case h{)Q ee quark & l instead of 1 — l/4({f — v') 2 . Wigner rotation 
contributions cancel also each other, which is seen to be due to the unitarity of Wigner 
rotations. 

Only the 1/4 coming for the active quark current has no counterpart in the r 2 's. One 
notes once more that this 1/4 is present in the j° matrix element of any elastic transition. 
It will be present also for a free quark. Let us also recall that the absence of 1/4 in baryon 
sum rules has nothing physical : it is due to the definition of invariant form factors, which 
are differently related to j° matrix elements. In terms of the latter, the only difference 
between mesons and baryons comes from the fact that % sums run over two spectator 
quarks instead of one ; but this refers to what we can term the compositeness contribution 

P space P wigner' 

In the following two subsections, because of the BIW sum rule, one needs only 
discuss p 2 , from which the parallel comments on the inelastic contributions can be deduced 
trivially. 



3.1 Non relativistic expansion 

It is useful to comment briefly on the order of magnitude of the contributions in a non- 
relativistic expansion, i.e. in powers of v/c where v/c is now the internal velocity of the 
light quarks. We know that this velocity is not actually small ; nevertheless, this expansion 
gives more physical insight. The dominant contribution to p 2 is the first one p 2 pace , from 
the spatial wave function. One finds at lowest order : 

P 2 « P%acc « m 2 (0|^ 2 |0) = O {{v 2 /c 2 )- 1 ) (41) 

_ m 2 R 2 - m ^ Q k Q case This corresponds for r to the use of the usual dipole formula. 
The Wigner rotation contribution is on the contrary highly suppressed : 

0^ = 0(v 2 /c 2 ) . (42) 

Pdirac = 1/4 i s * n between : O ((t> 2 /c 2 ) ) ; then it is actually the dominant relativistic effect 
coming from spin. This contrasts with Ref. ; there, the discussion identifies the spin 
effect with the Wigner rotation effect only. Denoting ground state averages as ( ) , we can 
also write a v/c expansion of the full p 2 : 



2 1 I P 2 

Hungrier 24 I m 2 
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P 2 = ™V)o + \ (jf V + z 2 p 2 ) Q + \ + 0(vVc*) (43) 

where the second term is 0((t> 2 /c 2 ) ) and comes for the "spatial" contribution p 2 space . It 
is seen that the lowest order relativistic corrections cannot be reduced in general, to the 
famous 1/4, as seems often assumed. Moreover they do not come from spin only ; 

Place = m 2 (z 2 ) o + 1/2 (p V + z 2 p 2 ) o + \ + 0{vV^) (44) 

contributes to them as observed in Ref. ||. Finally these additional relativistic corrections 
are not independent of the potential. Indeed (p 2 z 2 + z 2 p 2 )o has not even a definite sign. 

Truly, it is a convention to take m 2 (z 2 )o as the starting point of the v/c expansion. 
Indeed, it must not be forgotten that the average are taken on internal wave functions, 
which in realistic cases shall include by themselves relativistic binding corrections. Conse- 
quently, since eq.(|4"3"D does not make explicit the latter relativistic effects, it is not the full 
v/c expansion around a truly nonrelativistic limiting case. Nevertheless, the relativistic 
corrections written in eq. (|43"D may be considered as a reasonable answer to the question 
inasmuch as we are concerned only with relativistic corrections to p 2 due to the center-of- 
mass motion of hadrons. 



3.2 The lower bounds 



It is obvious that in this model there is no upper bound since p 2 can be arbitrarily large 
in the nonrelativistic regime, e.g. m 2 fl for an harmonic oscillator. On the other hand, we 
have found in the preceding paper [1[] a lower bound Inf p 2 = 3/4, which is larger than 1/4 , 
the famous Bjorken lower bound. To investigate further the compositeness contribution 
P%ace + Pwignev ft reveals useful to discuss the meaning of this bound. One may wonder 
why one cannot reach 1/4. Indeed, it would be tempting to suppose, by reverting the above 
argument, that one could reduce arbitrarily the "compositeness" contribution by reducing 
m 2 (z 2 )o, which corresponds to going to a highly relativistic situation (R 2 — > 0). But the 
idea is wrong, due to the relativistic effects. Part of the reason is the Wigner rotation 



contribution p 



wigner i 



which, for \p\ ^> m, amounts to 1/6. But another reason lies in the 



corrections to the "spatial" contribution p 



space ' 



Indeed, when m 2 (z 2 )o is reduced by going to large average momenta, p° is increasing 



in average and then, the final outcome is that the full expression of p 



2 

space ' 



( |(- 



-zp 



is bounded from below. We can see this by the change of variable introduced in 
as dx = p -^ : 





such 



p/m = sh x , p°/m = ch x 



(45) 



x 



Argsh p/m = Argch p°/m 



x coincides nonrelativistically with p, but behaves as logp/m for large p 

Write Plpace aS : 



(46) 

\p\. We can 



p°z + zp° 



(47) 
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where <po is the space part of the wave function. Then, using rotation invariance of <po 



p°z + zp° 



-<Po 



( d 1 p \ 
v dp 2 p / 



Tx + \ th x ) «"> 



Since d 3 p = An sh 2 x ch x dx on spherically symmetric wave functions, we define 



(48) 



g = v^4~7r sh xV ch x<£> (m sh x) 



(49) 



one then finds 



' d 1 p\ 

Po~r + «— 
v dp 2p ) 



dx 



thx) g 



J dx g ^1 



dx 7 



9=1 + 



dg 



dx 



(50) 



The second term is positive and has as lower bound, (obtained by spreading indefinitely 
g(x)). Therefore 



Inf 



p z + zp 







1 

3 



(51) 



Note that <po can be chosen arbitrary because the mass operator is arbitrary, except that 
Mq > 0. Since the Wigner rotation contribution is positive and smaller than 1/6 (since 
p/Po + rn <1), one gets : 7/12 < Infp 2 < 9/12 = 3/4. The more complete argument of |l] 
is necessary to show that : 

Inf p 2 = 2 (52) 

exactly, i.e. the Wigner rotation attains its maximum at the lower bound of p 2 . This 
reflects the fact that the bound is attained at large average of \p\. 

Since p 2 space is seen to be a major contribution to the lower bound 3/4 = 1/3 + 1/6 + 
1/4, it is worth specifying its origin. First it is an effect of compositeness, along with Wigner 
rotations. More specifically, it reflects the pure effect of spatial extension of the bound state. 
Finally, still more specifically, the bound 1/3 in eq. (|5T| ) reflects the effect of the relativistic 
transformation on the spatial wave functions. Indeed, nonrelativistically p 



space 



(m 2 (z 2 ) ) 

could be made arbitrarily small . The effect of the relativistic transformation law amounts 
practically to replacing m by the larger p° (as noticed in 0]). Then the contribution is 
bounded from below as just demonstrated, and in accord with the qualitative guess p° ~ |p| 
for large |p| and |p||-z| ~ 1. 



4 Conclusion 

In conclusion, we see that the heavy quark limit of quark models for currents based on 
the Bakamjian-Thomas construction of states exactly satisfies the important sum rule or 
duality relation discovered by Bjorken and further analyzed by Isgur and Wise. This 
property is essentially due to the nontrivial fact that the wave functions in motion satisfy a 
closure relation. This, added to the previous demonstration that it is covariant and scales 
as required by the Isgur- Wise heavy quark symmetry relations, increases the interest of the 
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model. Another aspect illustrated by the present analysis is the capacity of quark models 
to give a physical insight in the saturation of general field theoretic relations through the 
use of the concepts of bound state physics. Finally, one is stimulated to investigate the 
subdominant regime in 1/rriQ which is involved in other important sum rules such as the 
Voloshin "optical" sum rule ; in general, properties such as covariance or the conservation 
of the current are lost, but could perhaps be recovered under certain conditions. 
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